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ABSTRACT 

We demonstrate that the comparison of Tully-Fisher relations (TFRs) derived from global 
H 1 line widths to TFRs derived from the circular velocity profiles of dynamical models (or 
stellar kinematic observations corrected for asymmetric drift) is vulnerable to systematic and 
uncertain biases introduced by the different measures of rotation used. We therefore argue that 
to constrain the relative locations of the TFRs of spiral and SO galaxies, the same tracer and 
measure must be used for both samples. Using detailed near-infrared imaging and the circular 
velocities of axisymmetric Jeans models of 14 nearby edge-on Sa-Sb spirals and 14 nearby 
edge-on SOs drawn from a range of environments, we find that SOs lie on a TFR with the 
same slope as the spirals, but are on average 0.53 ±0.15 mag fainter at /Ts-band at a given 
rotational velocity. This is a significantly smaller offset than that measured in earlier studies of 
the SO TFR, which we attribute to our elimination of the bias associated with using different 
rotation measures and our use of earlier type spirals as a reference. Since our measurement 
of the offset avoids systematic biases, it should be preferred to previous estimates. A spiral 
stellar population in which star formation is truncated would take « 1 Gyr to fade by 0.53 mag 
at A'5-band. If SOs are the products of a simple truncation of star formation in spirals, then 
this finding is difficult to reconcile with the observed evolution of the spiral/SO fraction with 
redshift. Recent star formation could explain the observed lack of fading in SOs, but the offset 
of the SO TFR persists as a function of both stellar and dynamical mass. We show that the offset 
of the SO TFR could therefore be explained by a systematic difference between the total mass 
distributions of SOs and spirals, in the sense that SOs need to be smaller or more concentrated 
than spirals. 

Key words: galaxies: elliptical and lenticular, cD — galaxies: evolution — galaxies: kine- 
matics and dynamics — galaxies: spiral — galaxies: structure 



1 INTRODUCTION 

Both early- and late-type galaxies follow scaling relations. In the 
case of elliptical galaxies, the radii, velocity dispersions and lu- 
minosities lie on a Fundamental Plane ( [Djorgovski & Davis|1987[ 
[Dressier et al. ll987i. The luminosities of disc galaxies scale with 
their sizes JFreeman 197 0), rotation speeds ( |Tully & Fisher|1977| ) 
and colours ( Tully et al. 1982 ). The existence, scatter and evolution 



of these scaling relations provides powerful constraints for models 
of galaxy formation and evolution. 

In this work we concentrate on the Tully-Fisher relation (TFR), 
first discovered by [Tully & Fisher] ( | 1977) 1. In its original form, the 
TFR relates the global H I line widths of disc galaxies to their total 
luminosities. The line width of a rotationally-supported galaxy is a 
measure of the difference between the maximum rotation velocities 
of the approaching and receding sides. A galaxy's luminosity is re- 
lated to the mass of its luminous component (and the evolutionary 
state of the stellar populations). The observed correlation between 
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global line widths and luminosities is therefore a manifestation of 
a more general connection between galaxy masses and rotational 
velocities. A correlation is of course expected for gravitationally 
bound systems in which the mass distribution is dominated by lu- 
minous matter. 

Empirically, the correlation is tight enough to allow the TFR to 
be inverted and used to determine distances (e.g. [Sakai et al.|2000[ 
[Tully & Pierce|2000[ hereafter TPOO), but systematic variations of 
the residuals can arguably be used to infer the properties of dark 
haloes l [Courteau & Rix|1999^ . Efforts to simultaneously reproduce 
the zero-point of the TFR and the galaxy luminosity function with 
semi-analytical models of galaxy formation may reveal further in- 
formation on the role of the halo, the importance of feedback and 
the evolution of the stellar populations (e.g. [Baugh[[2006[ [Duttonj 
[et al.|2007| l. It is therefore clear that investigating variations of the 
local TFR with galaxy luminosity, size, type and other parameters 
should provide valuable insight into the growth of luminous struc- 
tures in the Universe. 

There are systematic variations in the slope, intercept and scat- 
ter of the TFR as a function of galaxy type. [Roberts] (fl978j and 
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[Rubin et al.H1985^ found that Sa spirals are fainter than Scs for 
a given rotational velocity at optical wavelengths. At least some 
of this difference has been attributed to variation in the shape of 
rotation curves related to the bulge-to-disc ratio (Verheijen 200 1[ 
|Noordermeer & \^ rheijen 2007]l. The difference is much smaller 
in the near-infrared ( .Aaronson & Mould] 1983[ [Peletier & Willner] 
[I993J. This has been taken as evidence that the discrepancy at opti- 
cal wavelengths is at least partially due to the effects of recent star 
formation on the luminosity, which is reduced in the near-infrared. 
This effect has been measured for a large, bias-corrected sample 
by [Masters et al.| ( |2008[ l. They find no evidence of any difference 
between the zero points of the Sb and Sc TFRs and a small effect 
for Sa galaxies, which they find have a variable offset compared to 
Sb-Sc spirals: ~ 0.2 mag brighter at v~ 150 km s^', no difference 
at V ~ 180 kms^' and ~ 0.4 mag fainter at v ~ 250 kms^'. 

Particular attention has been paid to SO galaxies, which are the 
earliest galaxies with rotationally-supported stellar discs. The fact 
that they are more common in the centres of clusters, rarer on the 
outskirts, and rarest of all in the field, the opposite of the trend ob- 
served for spirals, suggests that environmental processes transform 
spirals into SOs ( [Spitzer & Baade[|T95T| [Dressler| fT980'. 'Dressier 
etal.[1997). Processes such as ram-pressure stripping ( Gunn & Gott 



1972) , strangulation ( [Larson et al.f l980i and harassment ( jMoore 
et al.[[T996) may have removed the gas from the disks of spiral 



galaxies and left them unable to form stars. SOs do however exist 
in the field, where such processes should be insignificant. Passive 
evolution or the effects of active galactic nuclei in spirals have been 
suggested as possible formation mechanisms outside dense environ- 
ments ( [van den Bergh[2009^ . Problems remain with the disk fading 
picture, however: the bar fractions of spirals and SOs are discrepant 
([Aguerri et al.|2009{ [Laurikainen et al.|2009), in a given environ- 



ment SOs are at least as bright as spirals i Burstein et al.|2005) , and 
the bulge-to-disc ratios of SOs are difficult to reconcile with those 
of their presumed progenitors jChristlein & Zabludoff|2004) . See 
[Blanton & Moustakas[j2009) for a modem review of this issue. 

Unless they are violent, these processes should make only a 
slight change to the kinematics and dynamical masses of SOs, but 
a significant change to their luminosities. Assuming they do not 
have systematically different dark matter fractions, this change in 
stellar population should cause SOs to have fainter luminosities for 
a given rotational velocity and lie offset from the TFR of spirals. 
The size of any offset and the difference between the scatters of the 
SO and spiral TFRs can therefore be used to constrain models of SO 
formation and evolution ( [Bedregal et al.|2006! hereafter BAM06). 
A well-constrained SO TFR also raises the possibility of extending 
the TuUy-Fisher distance determination technique to earlier galaxy 
types, and of probing variations of mass-to-light ratio with redshift. 

Measuring the TFR of SOs is difficult because they do not usu- 
ally have extended gaseous discs. Instead of using a gaseous emis- 
sion line width as a measure of maximum rotation, several groups 
have used resolved rotation curves derived from stellar absorption 
lines. These raw observations must then be corrected for asymmet- 
ric drift (i.e. pressure support) before they can be presumed to trace 
the circular velocity (and therefore the potential of the galaxy). A 
number of groups have done just this and found that SOs are indeed 
fainter than spirals for a given rotational velocity. [Neistein et al.| 
( [199 9^ found that SOs lie on a TFR offset by 0.5 mag at /-band from 
the spiral TFR (with 0.7 mag intrinsic scatter). [Hinz et ar| ( [2001[ ) 
and |Hinz et al.[p"003| l found an offset of 0.2 mag at / and //-band 
(with 0.5-1.0 mag intrinsic scatter). [BAM06[ combined these data 
with new observations of galaxies in the Fornax cluster and found 
that the combined sample of 60 SOs lies around 1 .5 mag below the 



spiral TFR at S-band (0.9 mag scatter) and 1.2 mag below the spi- 
ral TFR at ATj-band (1.0 mag scatter). A promising alternative to 
stellar absorption line kinematics is sparse kinematic data from the 
emission lines of planetary nebulae. A pilot study of the nearby SO 
galaxy NGC 1023 with the Planetary Nebula Spectrograph found 
that this particular SO lies around 0.6 mag below a spiral TFR at 
A:5-band jNoordermeer et al.|2008l ). 

An alternative approach to measuring and correcting the ro- 
tational velocities of the earliest type galaxies is to construct dy- 
namical models that, in principle, directly probe the circular veloc- 
ity of a system. At the expense of model-dependent assumptions, 
this approach removes the uncertainties associated with the asym- 
metric drift correction, which may be both significant and system- 
atic. |Mathieu|erariJ2002l did this for a sample of six edge-on SOs 
and found that they lie 1.8 mag below an /-band spiral TFR (with 
0.3 mag intrinsic scatter). Modelling also allows the TFR to be ex- 
tended to elliptical galaxies with little or no rotational support. As 
with SOs, previous dynamical modelling studies of ellipticals have 
found that they are fainter than spirals for a given rotational ve- 
locity: [Franx[ (1993), 0.7 mag at /?-band; [Magorrian & Ballantyne| 
( [2OOT] ), 0.8 mag at /-band; [Gerh ard et al.'f2001!l, 0.6 mag at R- and 
1.0 mag at g-band; [De Rijcke e t al. (20071, 1.5 mag at B-, 1.2 mag 
at /f^-band. 

However, all previous studies comparing the TFRs of SOs and 
ellipticals to those of spirals have used independently determined 
reference TFRs for spiral galaxies (taken from, e.g., [Sakai et al.| 
2000, TPOO). These reference TFRs are derived in the traditional 
way, using global line widths from Hi observations. In Section |2] 
we discuss the technical and observational issues involved in char- 
acterizing the enclosed mass of a disc galaxy in an unbiased way. 
We argue that comparisons between TFRs derived from different 
measures of rotation are not a priori justified and are likely to intro- 
duce artificial and uncertain offsets in practice. 

The main astrophysical goal of the present work is thus to 
compare the TFRs of spiral and SO galaxies in a way that is not 
vulnerable to these uncertain systematics, by using measures of lu- 
minosities and rotational velocities determined identically for both 
SOs and spirals. Readers who are more interested in the results of 
this comparison and its implication for the evolution of SOs, may 
wish to skip ahead to Section [3] where we present our sample. In 
SectionHlwe describe the numerical method used to simultaneously 
determine the TFRs of spirals and SOs. We then present TFRs for 
the two samples in terms of near-infrared and optical luminosities, 
stellar mass and dynamical mass. In Section [5] we discuss our re- 
sults in the context of models of SO formation and evolution and 
previous analyses of SO TFRs. We conclude briefly in Section[6] 



2 MEASURES OF ROTATION 

2.1 The problem 

To construct a TuUy-Fisher relation one needs a measure of rota- 
tion. This can be derived from, for example, the global Hi line 
width, spatially-resolved H I, Ha or stellar rotation curves or veloc- 
ity fields, sparse tracers (e.g. globular clusters and planetary nebu- 
lae), or the circular velocity of dynamical models at some fiducial 
location. In principle, all these measures are related to the rotation 
of the galaxy, but the relationship of these single numbers to the en- 
closed mass is complicated by the fact that observed rotation curves 
(and indeed modelled circular velocity curves) do not always flat- 
ten. Moreover, it not obvious that these observational measures di- 
rectly and accurately comparable to each other. Both the nature of 
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the observations and the corrections applied for different techniques 
could conceivably introduce systematics that vary with galaxy prop- 
erties such as type, mass or size. 

For example, the transformation of a global H I line width to 
a measure of rotation involves corrections for instrumental broad- 
ening and turbulent motion (e.g. |TuUy & Fouque|1985[ |Verheijen| 
|2001^ . The extent to which these corrections are able to recover 
the intrinsic motion of the tracer or the true circular velocity of 
the galaxy potential has recently been questioned (Singhal 2 008^ . 
Meanwhile, stellar kinematics must be corrected for asymmetric 
drift, which typically involves making approximations which are 
difficult to justify prima facie (see Section |2.2^ . Finally, the con- 
struction of dynamical models often involves simplified modelling 
of the full distribution function (e.g. [Mathieu et al.|[2002^ or as- 
sumptions about intrinsic morphology, anisotropy and the role of 
dark matter (e.g. |Williams et al.|2009| hereafter WBC09). 

These considerations place the direct comparisons of TFRs de- 
rived using different measures of rotation on shaky ground. Offsets 
which are ascribed to fading or brightening could in fact be due to 
systematic biases that differ between the measures of rotation used. 
For typical slopes of the TFR, a relatively small systematic differ- 
ence in velocity of ~ 0.1 dex introduces an offset in the TFR that is 
indistinguishable from a --^ 1 mag difference in luminosity. This is 
of the order of the typical zero-point offsets claimed for SO TFRs 
( |Neistein et al.|1999|[Hmz et al.|200l|[2003lpathieu et al.|2002| 
|Bedregal et al.|2006 ll, so it is crucial that we eliminate (or at least 
quantify) the systematics introduced by such comparisons. In the 
remainder of this section, we examine the differences between the 
results of three different measures of rotation when applied to the 
same 28 early-type disc galaxies, and discuss the implications for 
comparisons of the zero-point of the TFR. 



2.2 Four possible measures of rotation defined 

We aim to compare measures of rotation derived from global H I 
line widths (which we denote vjji), resolved gas rotation curves 
(vgas), stellar kinematics (vjjiift) and dynamical models (vj^odel) for 
the 28 galaxies in the present work. This sample, which consists of 
14 spirals and 14 SOs, is described in SectionlS] 

For I'm, the rotational velocity derived from global Hi line 
widths, we adopt the quantity vrot from HYPERLEDA ( |Paturel| 
|et al.|2003| l. This is the inclination-corrected maximum gas rotation 
velocity, based on an average of independent H I line width determi- 
nations taken from the literature. Strictly speaking, HYPERLEDA 
also uses spatially-resolved Ha rotation curves to calculate vrot, 
but these observations were only present in HYPERLEDA for four 
of the galaxies in our sample, and there was no evidence that they 
systematically disagreed with the H I line width values at any more 
than the 0.01 dex level. 

For Vgas, we fit Gaussians to the position- velocity diagrams 
(PVDs) presented in |Bureau & FreemanH1999| l. These are derived 
from [N II] emission lines. We only use the region of the PVD where 
the rotation curve is flat. This restriction means that it is not possible 
to measure Vgas for many of the galaxies in our sample, either be- 
cause emission was not detected, or because it was not sufficiently 
extended. 

We calculate the asymmetric drift-corrected stellar velocity, 
Vdrift, using the observed stellar kinematics presented in [Chung & | 
[Bur eau (2004). In order to more accurately measure the local circu- 
lar velocity we attempt to correct these observations for the effects 
of line-of- sight integration and asymmetric drift using the same 
recipe as used in[BAM06] which we describe briefly below. We 



make no attempt to improve this recipe, since one of our goals is 
to assess the biases in previous determinations of the TFR based on 
stellar kinematics. The raw observations and their corrected values 
are shown in Fig.fTl 

The line-of-sight correction for an edge-on disc assumes that 
the local azimuthal velocity v^{R) at a galactocentric distance R 
(in the galaxy's cylindrical coordinate system {R,(j),z)) is related to 
^LOsi^)^ the luminosity- weighted mean line-of-sight velocity at a 
projected distance from the galaxy centre x', by the equation 



vlos [X , 



JoV^{R)p{R)cos(t>dz' 



(1) 



Following the notation of, e.g., [Binney & Tremaine| pOOSl, 
{x',y' ,z') is the observer's coordinate system on the sky, in which 
thex' axis is aligned with the projected major axis of the galaxy and 
the z' axis is along the line of sight, p (R) is the luminosity density. 
The integral is evaluated along that line-of-sight z' ■ By assuming 
that the local streaming velocity v^ is independent of R, one can 
use the above relationship to infer it. This assumption is valid only 
if/? is at least the radius at which the intrinsic circular velocity curve 
flattens. To evaluate these integrals we assume p{R) is an exponen- 
tial disc of scale length R^, which we measure from the major-axis 
surface brightness profiles of our Ks-band images at radii outside 
the bulge (Bureau et al.|2006| l. 

We then apply a correction for asymmetric drift to v^ (R) to 
derive v^ (inft(/?), an estimate of the true local circular velocity: 



"c, drift V 



■4w 



-oi(R) 



'Rd 



(2) 



where CT0 {R) is the azimuthal stellar velocity dispersion. This equa- 
tion is derived from the Jeans equations (Jeans'1 922| [Binney &] 
|Tremaine[2008} by assuming a thin disc in a steady state and the 
epicyclic approximation (ai/o^ = 0.5 for constant V0). To esti- 
mate aij,{R), we assume that it is equal to an exponential fit to 
the observed stellar line-of-sight velocity dispersion Ohosix')- This 
assumption is of course flawed, since the observed dispersion in- 
cludes a contribution from the changing component of the az- 
imuthal velocity along the line-of-sight. This effect will bias the 
asymmetric drift correction to be too high. However, the error is 
small in practice, and, crucially for the present work, it is not a 
strong function of galaxy type in the range SO-Sb. Under the as- 
sumptions described above, for disks with typical levels of pres- 
sure support (v^/<J^ ~ 5 at the relevant radii ~ 2-4R^), numerical 
integration shows that the error introduced by wrongly assuming 
0,j,{R) = Olos(j^') causes Vc^jift to overestimate the true circular 
velocity by ~ 0.04 ± 0.01 dex. This error is almost independent of 
v^/cr^ in the range 2 < v^/a,^ < 20, a quantity which is in any case 
independent of morphological type for our sample. 

Finally, the local circular velocity of the mass model, v^ model) 
is derived from the dynamical models presented in |WBC09] In that 
work, we modelled the mass distribution of each galaxy by assum- 
ing it is composed of an axisymmetric stellar component with a con- 
stant mass-to-light ratio and a spherical NEW dark halo ^|Navarro| 
|et al.|199Tj l. The axisymmetric stellar component is based on an an- 
alytic fit to the observed near- infrared photometry of the composite 
bulge and disk. The absence of a triaxial bar component in our mod- 
els does not significantly affect either the global properties of the 
models or the local kinematics at the large radii relevant to the TFR 
(see section 5.5.3 of |WBCO^ for a more extensive discussion of this 
issue). We assumed a particular relationship between halo mass and 
concentration jMaccio et al.|2008] >. We determined the three param- 
eters of the models (the stellar mass-to-light ratio, virial halo mass 
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Figure 1. Observed, corrected and modelled kinematic data for all sample galaxies. The black points are the observed mean line-of-sight stellar velocities vlqs 
and the grey points the line-of-sight velocity dispersions Olos. both from Chung & BureauH2004^ . The open red squares are estimates of the circular velocity 
derived by correcting the observations for line-of-sight and asymmetric drift effects (v^ dnft). There are missing points at small radii because the correction is 
unreliable where the local dispersion is large compared to the line-of-sight velocity. The solid blue line is the circular velocity of the galaxy dynamical model 
from|WBC09|(Ve model)- These models consist of contributions from a stellar component (dashed blue line) and a dark halo (dotted blue line). The vertical line 
is the radii beyond which the observed rotation curve is flat, which is important for our calculation of Vdnft and I'model. the single values used to characteiize 
rotation. These are determined as described in the text. Their values for each galaxy are shown as filled red and blue squares respectively, at arbitrary large 
radii. 
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and orbital anisotropy) by comparing tlie observed stellar second 
velocity moment Vims (x') = [v^qs (^ ) + ''los (''^)] along the ma- 
jor axis to that predicted by solving the Jeans equations assuming a 
constant anisotropy (parametrized as j8- = 1 — CT- /aj). The solution 
of the Jeans equations under these assumptions is described in |Cap-| 
pellari ( 2008b. Note that while Ve.drift described above is weakly af- 
fectedby our assumptions that ctI/ctJ = 0.5 and Oip (R) = CTlos (x'), 
Vc,model suffers from no such biases. Our solutions of the Jeans 
equation are constrained by a fit to the observed second velocity 
moment along the line-of-sight. The second moment is strictly in- 
dependent of oi/a^ so no assumption for this ratio is required, 
and the implementation rigorously accounts for line-of-sight inte- 
gration. 

To construct a TFR, the radial profiles of the asymmetric-drift 
corrected stellar kinematics, V(..drift and of the model circular ve- 
locity, Vc^odel must be characterized by a single number, Vdrift and 
Vfjiodel respectively. We do this by determining the flat region of the 
observed rotation curve by eye. We then take Vinii to be the mean of 
the values of Ve_drift beyond this point. The continuous Jeans mod- 
elled circular velocity curve is evaluated at the same radii as the ob- 
served stellar kinematics. Its mean in the flat region of the observed 
rotation curve gives VmodelOf course, it is possible to evaluate the 
circular velocity of the models at arbitrary radii, but by restricting 
ourselves to those radii for which we have stellar kinematic obser- 
vations, we avoid using the mass models in regions where they are 
not constrained by the data. Fig. [T] shows these calculations for all 
sample galaxies. 



2.3 Comparison 

A comparison of the four measures of rotation is shown in Fig. [2] 
The upper panel demonstrates that (log(v^o(iel/^'drift)) = — 0.015 ± 
0.003 dex with an rms scatter of 0.019 dex for both spirals and SOs. 
This small difference, which is on average less than 4 per cent, is 
probably due to the additional approximations and assumptions in- 
troduced in the asymmetric drift correction compared to the Jeans 
models. We note in particular that it is consistent with the error due 
to assuming (T^(x') = cTlos(-'^)! as discussed in Section 2.2 The 
absolute offset between these measures is not significant in the con- 
text of the TFR and, crucially, there is no evidence of any systematic 
difference between SOs and spirals. The scatter is also extremely 
small. Since both Vdrift and Vmodel are derived from the observed 
stellar kinematics, their close agreement is not of itself independent 
proof of the accuracy with which they trace the true circular veloc- 
ity. The agreement of Vdrift with the results of detailed and robust 
Jeans modelling does suggest, however, that the approximations in- 
herent in eqns. (1-2) yield an approximately valid solution of the 
Jeans equations. 

The middle panel demonstrates fairly good agreement between 
the two measures based on gas, v^i and Vg^s- (log(vHi/vgas)) = 
0.032 ± 0.009 dex with a rms scatter of 0.029 dex. There is some 
evidence that vhi systematically exceeds Vgas by 5-10 per cent. 
Statistics are poorer because gas is rare in SOs by definition, but 
there is no evidence of any systematic difference between I'm/vgas 
in SOs and spirals. 

However, the bottom panel of Fig. [2] shows that v,j,odel (and 
therefore Vdrift) does not in general agree vm (and thus Vgas): 
(log(vHi/vnjodel)) = —0.094 ± 0.031 dex with a rms scatter of 
0.1 15 dex (or -0. 18 ± 0.043 dex with a rms scatter of 0.243 dex if 
the |BAM06l galaxies are included). Some of this difference and scat- 
ter is caused by pathological cases of I'm derived from single dish 
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Figure 2. Comparisons of the circular velocities of the mass models (v^oAa)^ 
observed stellar kinematics corrected for asymmetric drift (I'drift) and rota- 
tional velocities derived from global Hi line widths (vm) and ionized gas 
position-velocity diagrams (vVas). The solid blue triangles are the Sa-Sb 
galaxies from our sample and the solid red circles the SOs. The open red 
circles are the SOs presented in |BAM0'6| They do not estimate I'modcl so, in 
order to place these galaxies on the plot, we assume vjrif, « v^^jujei for them. 



measurements, i.e. spatially unresolved measurements of 21 cm sig- 
nal from galaxies with a warped gas disc, polar rings or gas-rich 
nearby companions that fall within the telescope beam (e.g. this is 
almost certainly the case for the SO galaxy NGC 1596, see [Chung] 
let al.|2006^ . 

But even when possible pathological cases are avoided, 
for example by ignoring SO galaxies, which are more likely 
to be problematic, there remains a clear systematic offset: 
(log(vHi/vmodel)) = -0.063 ± 0.005 dex (15 per cent) with 
0.02 dex rms scatter for the spirals in our sample. For these galax- 
ies, the rotational velocity derived from the global H I line width is 
significantly and systematically smaller than that derived from stel- 
lar kinematics or dynamical modelling. Because v^i is derived from 
spatially unresolved data, this could be explained by the global H I 
line width probing the local gas velocity out to different radii than 
our typical stellar kinematic data jNoordermeer et al.|2007^ . 

We find a similar systematic difference between v^jodel (or 
Vdrift) and Vgas: (log(vgas/vinodel)) = -0.086 ± 0.008. Vgas is de- 
rived from a spatially resolved PVD taken from the same long- 
slit optical spectrum with which we derive the stellar kinematics. 
We can therefore verify directly that it does not suffer from either 
pathological problems related to gas-rich companions, warps or po- 
lar rings, or doubt about the radii probed compared to our stellar 
kinematics. We show this in Fig.lSlfor the five sample galaxies with 
optical emission that extends to the disk. It is clear that the stars and 
gas are associated with the same galaxy and are observed at similar 
radii. Therefore, unless light-emitting gas is absent from the tangent 
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point along the line-of-sight in these edge-on galaxies (a possibility 
we cannot rule out), the emission line gas kinematics and absorption 
line stellar kinematics should imply the same local circular velocity 
and the same single characteristic measure of rotation, (vgas, Vjjjjft 
or V^odel)- 

But in these five galaxies galaxies, the stars and gas rotate 
with essentially the same velocity as a function of position in the 
outer, flat parts of the observed rotation curve. The stars have non- 
negligible dispersion so we know they do not trace the local cir- 
cular velocity by an amount that can be calculated using the Jeans 
equations. Given this observation, it is not surprising that Vgas is 
significantly less than v^jodel- 

In summary, while agreement between Vga., and vjji is fairly 
good, and from this we conclude that Vmojel and v^rift (based on 
stars seen in absorption) do not imply the same circular velocity as 
Vhi and Vg^s (based on gas seen in emission). This is not a new re- 
sult. In early-type disc galaxies, the observed stellar kinematics and 
the well-understood Jeans equations often imply a local circular ve- 
locity greater than the observed rotational velocity of gas (e.g. |Kent| 
1988; Kormendy & Westpfahl 1989; Bertola et al. 1995; Cinzano 
et al..l999,,Corsini et al.„1999,,Cretton et al..2000,,Vega Beltran, 



et al.|200l||Pizzella et al.|2004||Krajnovic et al.||2005[|Weijmans 



et al.|2008||Young et al.|2008^ . This puzzle is made more acute by 
the small velocity dispersions observed in the emission line PVDs. 
We will return to the gas dispersion in a future work. 

A definitive resolution of this issue is outside the scope of the 
present work and, fortunately, it is unnecessary for our goals. For 
the purposes of comparing TFRs (or for determining relative dis- 
tances; see e.g. jCourtois et al.|2009^ , it does not matter which mea- 
sure is 'right', as long as we use the same measure for the two sam- 
ples and the measure chosen does not introduce a systematic bias 
between samples. We see no evidence of this in our data, for any of 
the four measures explored. 

Our conclusions in this section have three consequences for 
previous works that measure the TFR of SO galaxies using stellar 
kinematics corrected for asymmetric drift (e.g. |Neistein et al.|1999[ 
|Hinz et al.|200T|[2003||BAM06[ >: (i) the asymmetric drift correc- 
tion they use yields similar results to detailed Jeans modelling; (ii) 
the drift correction does not seem to introduce a systematic bias be- 
tween spirals and SOs and, if applied to samples of both spirals and 
SOs, can indeed be used to compare the TFRs of the two classes; 
(iii) however, a TFR derived from asymmetric drift-corrected stel- 
lar kinematics cannot be directly compared to TFRs derived from 
global H I line widths or resolved emission line PVDs, as these pre- 
vious authors did. As we argue above, such comparisons are not a 
priori justified and, in the case of both our sample and the larger 
sample of [BAM06] they appear to introduce systematic errors that, 
on a Tully-Fisher plot, are of the same order (and in the same di- 
rection) as the TFR offsets that are usually interpreted as luminos- 
ity evolution. The conclusions of previous works regarding the off- 
set of the SO TFR should therefore be treated with caution. In ad- 
dition, given the large scatter we observe in log(vHi/vn,o(jel) and 
log(vHi/vcirift)' it is far from clear that any simple additive or mul- 
tiplicative correction could eliminate the bias. 

For our sample, we are able to side-step this issue completely 
by using a single measure of rotation for both spirals and SOs. In this 
way we compare the TFRs of SO, Sa and Sb galaxies in a manner 
which is not subject to the possible systematics discussed above. We 
adopt Vmodei as our fiducial measure of rotation for both spirals and 
SOs. The figures and discussion that follow use only this measure. 




-400 L 



200 



Figure 3. Compaiisons of [Nil] A6584 position-velocity diagrams, stel- 
lar kinematics, and the circular velocity implied by Jeans modelling for 
the five galaxies in the sample with optical emission that extends to the 
disc. From top to bottom: ESO 240-Gll, IC 5096, IC 5176, NGC 5746 and 
NGC 6722. The filled contours are the continuum-subtracted optical spec- 
tra originally presented in [Bur eau & Freeman 1 19991. The observed mean 
stellar velocities jChung & BureaU|2 004i and the circular velocity profiles 
inferred from detailed Jeans modelling of the light distribution and stellar 
kinematics jWBC09f are overplotted as red crosses and blue lines respec- 
tively. 



3 SAMPLE AND DATA 

To make the reliable comparison we described in the previous sec- 
tion, we use the sample of 14 early-type spirals (Sa and Sb) and 
14 SOs presented in Table [T Dynamical models of these galaxies 
were constructed in|WBC09 All of the galaxies are oriented close 
to edge-on (within ~ 5°). Many of the galaxies have boxy bulges, 
which are believed to be bars viewed side-on (e.g.| Kuijken & Mer-| 
rifield 1995; Merrifield & Kuijken 1999. Bureau & Freeman 1999) 
iChung & Bureau,2004p , and it was to probe these bulges that the 
sample was constructed and observed. The bars should not, how- 
ever, affect the results we present here because we use measures of 



© 2010 RAS, MNRAS OQO.fTpT] 



The Tully-Fisher relations of early-type discs 7 



Table 1. Sample. 



Galaxy 


Mk, 


Mb 


log(M*/M,.,) 


log(Mdyn/Mv,) 


Ri 


^'model 


Vdrift 


VHl 


^'gas 


2MRS LDC Clu 




(mag) 


(mag) 






(arcsec) 


(kms-') 


(kms-') 


(kms"') 


(kms-') 


group .size 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(V) 


(8) 


(9) 


(10) 


(11) (12 



NGC 128 


-25.35 


-21.40 


11.47 


ESO 151-G004 


-24.85 


-20.41 


11.40 


NGC 1032 


-24.47 


-20.71 


11.13 


NGC 1381 


-22.76 


-18.77 


10.49 


NGC 1596 


-22.86 


-18.98 


10.56 


NGC 2310 


-22.48 


-19.43 


10.25 


ESO311-G012 


-23.06 


-19.92 


10.45 


NGC 3203 


-23.89 


-19.99 


10.76 


NGC 3957 


-23.29 


-19.20 


10.76 


NGC 4469 


-23.01 


-19.02 


10.56 


NGC 4710 


-23.47 


-19.40 


10.56 


NGC 5084 


-24.73 


-21.85 


11.16 


NGC 6771 


-25.08 


-20.81 


11.40 


ESO 597-G036 


-25.36 


-20.70 


11.58 


NGC 1886 


-22.90 


-19.73 


10.47 


NGC 3390 


-24.90 


-21.49 


11.19 


NGC 4703 


-25.43 


-20.92 


11.38 


PGC 44931 


-24.46 


-21.33 


11.07 


ESO 443-G042 


-23.49 


-20.21 


10.76 


NGC 5746 


-25.64 


-22.30 


11.63 


IC 4767 


-23.69 


-19.72 


10.93 


NGC 6722 


-25.38 


-21.97 


11.49 


ESO 185-G053 


-24.22 


-20.09 


11.01 


IC 4937 


-24.67 


-20.46 


11.18 


NGC 7123 


-25.29 


-20.81 


11.35 


IC 5096 


-24.93 


-21.21 


11.27 


IC5176 


-23.31 


-19.95 


10.60 


ESO 240-GOl 1 


-24.79 


-21.56 


11.56 



11.56 
11.46 
11.22 
10.59 
10.66 
10.48 
10.53 
10.94 
10.85 
10.68 
10.68 
11.22 
11.46 
11.73 

10.63 
11.27 
11.49 
11.27 
10.90 
11.75 
10.97 
11.60 
11.06 
11.24 
11.40 
11.36 
10.72 
11.57 



SO 










16.7 


369 


379 






8.2 


301 


311 






33.0 


281 




284 


254 


12.3 


208 


220 






20.5 


247 


249 


98 




25.2 


171 


173 






28.3 


196 


192 


96 




16.3 


222 


228 






19.1 


213 


213 




178 


16.5 


180 


199 






23.6 


187 


199 


160 


140 


34.7 


303 


306 


310 




16.2 


337 


337 




284 


14.4 


345 


365 






Sa-Sb 










14.5 


178 


192 


155 




20.4 


257 


279 


210 


229 


29.2 


269 


272 


246 


229 


21.1 


242 


225 


200 


180 


11.4 


190 




166 




69.1 


365 


393 


319 


295 


8.2 


195 


219 






12.5 


323 


322 




263 


5.0 


244 


267 






18.7 


208 


210 






19.2 


327 


335 






14.5 


290 


300 


262 


244 


20.9 


194 


217 


164 


150 


35.7 


305 


306 


268 


243 



43 
30 



210 
13 
300 
300 
15 
97 
3 



210 
80 
80 
38 

45 
97 
97 
42 
3 

5 

3 



Fornax 



Hydra 

Virgo 
Virgo 



Hydra 



Notes. Sample divided into SOs and spirals based on classifications of |Jarvis|jl986) , rde Souza & Dos Anjos|{T987t , |ShawH1987| and [Karachentsev et al.| 
(1993). (1) Galaxy name. (2)-(3) Total Ks-hand and B-band absolute magnitude. (4)-(5) Total stellar and dynamical mass. (6) Disk scale length measured using 
Xj-band images presented in |Bureau et al.|{2006t . (7) Circular velocity derived from the best-fit mass model presented in |WBC09] (8) Rotational velocity 
derived from the observed stellar kinematics presented inl Chung & Bureau] j2004l, corrected for line-of-sight effects and asymmetric drift. Missing values are 
cases where the drift correction was too large to be considered reliable. (9) Rotational velocity derived from the global H I line width, i.e. the quantity vrot 
from HYPERLEDA where available. (10) Ionized gas rotational velocity measured by a Gaussian fit to [Nil] position-velocity diagrams presented in B ureauj 
|& FreemanH1999y (11) Number of companion group members in the Low Density Contrast 2MRS Group Catalog ( Crook etak2007i. Ellipses indicate that 
a galaxy was not found in the catalog, implying that it had 2 or fewer nearby companions. (") ESO 151-G004 is not present in the 2MASS Extended Source 
Catalog (Jarrett et al.|2000) and so is not included in the 2MRS database. (12) Cluster membership. For full details of how the quantities in Columns (2)-(10) 
were determined, see Section [T2] and Section[3] 



rotational velocity well outside the boxy bulge region. The fraction 
of galaxies with boxy bulges in the sample (~ 75 per cent) is not 
in any case significantly different to the fraction of barred galaxies 
in the local universe (~ 65 per cent; see, e.g. |Eskri"dge et al. 2000[ 
|Whyte et al.|2002|pm"nova & Jogee|2007| ). 

The sample was selected without regard for environment, and 
so includes galaxies in a wide range of environments. Two SOs are 
members of the Virgo cluster, an SO and a spiral are members of the 
Hydra cluster, and an SO is a member of the Fornax cluster. Cross- 
referencing our sample with the 2MASS Redshift Survey (2MRS) 
Group Catalog (Crook et al.|20 07 1 demonstrates that the remainder 
of the sample is made up of members of intermediate-size groups 
and relatively isolated galaxies. 

All the galaxies are relatively bright and fast-rotating, and the 
spirals and SOs cover the same range of luminosities. This means 
that we are not vulnerable to variation in the slope of the TFR at 
the low or high mass end, as has been observed (e.g. |Peletier &| 



|Willner|199"3||Noordermeer & Verheijen|2007| ). We further note that 
the sample contains some very bright SOs, such as ESO 151-G004 
and NGC 6771, that have large boxy bulges. If these bulges are in 
fact bars, and therefore the products of secular evolution, then these 
particular bright SOs are unlikely to be the elliptical-like products 
of major mergers, as has previously been suggested ( Poggianti et al.| 



[200T|p;fehrert et al.|2003[|Bedregal et al.|2008[[Barway et aL|2009) . 

In Fig. PI we present the spiral and SO TFRs as functions of 
A'^-band luminosity, B-band luminosity, stellar mass and dynam- 
ical mass, all discussed below. Throughout this work we assume 
an absolute solar magnitude Mks,q = 3.29 mag (Blanton & Roweis 
[2007} . 

Firstly, all the vertical axes in Fig. p] require a distance es- 
timate for each galaxy. As described in WBC09, we use surface 
brightness fluctuation (SBF) estimates where available (NGC 1381 
from [Jensen "erar][2003] and NGC 1596 from [Tonry erar][200T] ) 
and the Virgo cluster distance where appropriate (NGC 4469 and 
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Figure 4. Tully-Fisher relations for our samples of SO and spiral galaxies. 
From top to bottom, as a function of Vn,odel '■ total absolute Sj-band luminos- 
ity, total absolute B-band luminosity, total stellar mass and total dynamical 
mass. Spiral galaxies are shown as blue triangles and SO galaxies as red cir- 
cles. A median error bar is shown in the top left corner of each plot. The 
best fit relations found by inverse regression are shown as solid lines (blue 
for spirals, red for SOs). Their parameters are given in Tablel2] The dashed, 
dotted and dot-dashed lines are a selection of TFRs found by previous au- 
thors for different morphological types. See text for a full discussion of these 
relations. 



NGC 4710, |Mei et al.||2007| l. For all other galaxies we adopt 
the redshift distance presented in the NASA/IPAC Extragalactic 
Database, conected for a Virgocentric flow model ( [Mould et al.| 
j2000l. We assume a Wilkinson Microwave Anisotr opy Probe five- 



Komatsu et al. 



year (WMAP5) cosmology {Hq = 70 km s" ' Mpc^ ' 

|2009p . Of course, we avoid distance estimates derived from TFRs 

The total apparent magnitude at A'^-band is measured using 
the analytically derived total light in the multi-Gaussian expan- 
sions ([Emsellem et al.|1994[|Cappellari 



2002b of our near-infrared 



images, as described and tabulated in WBC09| The total appar- 
ent magnitude at S-band is taken from HYPERLEDA ( |Paturel| 
|et al.|[2003| l. We apply corrections for Galactic (Aq) and inter- 
nal (A;) extinction in both bands using the Galactic dust maps of 
[Schlegel et al. 1 1998|l and the internal extinction correction pre- 
sented in |Bottinelli et al.] ( fT995t , as tabulated in HYPERLEDA ^ 



|turel et aLl|2003^ , both of which we transform to the appropriate 
bands using the extinction law of |Cardelli et al.| ( |1989| ), implying 
th2i\.Ac{Ks)/Ac{B) =A;(Ks)/Ai{B) = 0.084. The internal extinc- 
tion correction is a function of both bulge-to-disc ratio and morpho- 
logical type. This correction is < 0. 1 mag at A'^-band for our spirals 
and smaller still for the SOs. 

With the exception of the Virgo galaxies and the two galaxies 
with SBF distances, we adopt characteristic distance errors corre- 
sponding to a 200 km s^' uncertainty in the flow model (e.g. Tonry 



|etal.|200^ , and 0.05 mag uncertainties in the apparent magnitudes. 
These uncertainties completely overwhelm any error in the Galac- 
tic or internal extinction at A^^-band. At S-band, we adopt uncer- 
tainties of 0.1 6 mag in Galactic extinction (Schle gel et al.|p~998[ l 
and 0.1 mag in the internal extinction. We emphasize, however, 
that internal extinction corrections in edge-on spiral galaxies are 
both large and uncertain. The i?-band TFR is problematic with even 
favourably inclined spirals, but with our edge-on sample the internal 
extinction correction and its uncertainties may introduce significant 
systematic errors. The B-band data are therefore included in this 
work merely to aid comparison with previous work, and we do not 
make use of these results in our interpretation. 

The stellar mass, M*, is derived from the total absolute 
ATj-band luminosity and the constant stellar mass-to-light ratio, 
{M/L)ks, of each galaxy determined in iWBC09 As described in 
Section [T2] this uses dynamical methods that make assumptions 
about the dark halo and stellar dynamics, but are not vulnerable to 
zero-point or colour-dependent uncertainties due to the initial mass 
function or stellar evolution. 

The dynamical mass, M^yn, is estimated using the total abso- 
lute A'^-band luminosity and the dynamical mass-to-light ratio at 
A'^-band for each galaxy, also determined in WBC09 in which it is 
referred to as {M/L)xg.noiMAo- Unless there is no dark matter (or its 
distribution closely follows that of the luminous matter), the dynam- 
ical (M/L) ratio is likely to be a quantity that varies significantly 
within a galaxy, increasing toward the edge of the optical disc where 
dark matter normally begins to dominate. The value we use is, how- 
ever, effectively an average for each galaxy within the radii con- 
strained by the stellar kinematic data, typically 2-3 Re(wheie R^is 
the projected A^^-band half-light radius), but dominated by the cen- 
tral regions. Because of this, Mjy,, is a somewhat ill-defined quan- 
tity. 

An alternative approach would be to adopt the mass within the 
three-dimensional half-light volume, Mj ;2- We show how this can 
be evaluated from observations or models of the circular velocity in 
disc galaxies in Appendix [A| We also show there that, in practice, 
2M 



\I2- 



'- Mlyn to within 13 per cent. Neither choice is optimal but 



in the context of the present work, we can think of no reason why 

© 2010 RAS, MNRAS OQO.fTpT] 



The Tully-Fisher relations of early-type discs 9 



one choice or the other would introduce systematic biases that vary 
between spirals and SOs. The choice is therefore somewhat arbitrary 
and has no affect on the discussion of the relationship between spi- 
rals and SOs that follows where we use M^y„. This follows precedent 
(it is the definition used by, e.g., Magorrian et al. 1998 Hiiri ng &| 
|Rix|2004] and |Cappellari et al.,2006| ) and avoids the possible con- 
cern that the argument we present in Section [F!2| is circular, but we 
emphasize that Mj m gives the same results and is a better defined 
quantity. 

For reference, we show in Fig. H] the spiral TFRs determined 
using global H I line widths at A'-band by TPOO, at A'^-band by|Mas-j 
Iters et al-lpOOS) ) and at S-band by TPOO and |Sakai et aTlpOOoV We 
do this assuming that the global H I line width, once corrected for 
inclination, broadening and turbulent motion, is exactly twice the 
rotational velocity. The zero points of the TPOO and |Sakai et al.| 
( |2000[ > TFRs were determined using known Cepheid distances to 
4 galaxies in the case of TPOO at A^-band (1 Sab, 2 Sb and 1 Sc; 
see the companion paper by |Rothberg et al.|2000^ , 24 galaxies for 
TPOO at B-band (3 Sab, 14 Sb-Sc and 7 Scd-Sd) and 21 galax- 
ies for |SakaietaT]pOOO) (3 Sa-Sab, 16 Sb-Sc and 2 Scd). We 
position the [Masters et al.| ( |2008| ) A's-band TFR on our plots using 
our adopted value of Hq. The difference between the two A^-band 
and two fi-band relations gives an idea of the uncertainty in the ab- 
solute location of the late-type spiral TFR (the difference between 
the A'-band filter used by TPOO and the 2MASS A'j-band filter is 
not significant in this context). These comparison spiral TFRs are 
unconstrained by observations above log(v/km s^') ~ 2.45, while 
a third of our sample is above this value. 

Unlike |BAM06| we do not shift the zero-points of the TPOO 
TFRs by -0.207 mag. |BAM06| did this to make the value of Hq im- 
plied by TPOO's data (77 ± 8 km s^ ' Mpc^ ' ) consistent with their 
adopted value (70kms^' Mpc^'). However, we believe that this 
correction was in error, since the absolute locations of the TPOO 
TFRs were set by independent Cepheid distances, which imply 
rather than assume a particular value of Hq. 

We also show an extrapolation to high mass systems of the 
baryonic TFR presented in |Stark et al.| ( |2009) , which was calibrated 
using gas-dominated dwarf galaxies (hence minimizing uncertain- 
ties associated with the stellar mass-to-light ratio). Neither of our 
mass estimates M* or Mjyn are designed to reproduce the bary- 
onic mass M\)^ used by [Stark et al.| ( (2009^ . In practice, however, 
one expects M* ~ M^ar to within a few per cent for gas-poor SOs 
and to within perhaps ten per cent for giant Sa-Sb spirals, while 
Mdyn ~ Mbai- to within a factor of 2 for both SOs and spirals (dark 
matter is no more than 50 per cent of the mass within the optical ra- 
dius of disk galaxies, e.g. |van Albada & Sancisi 1986, Persic et al. 
T996i[PaIunas & Williams|2000„Cappellari et al.|2006„Kassin eT^ 



locities. This effect would propagate into the rotation curves from 
which Vniodel and Vdrift are derived, and also into the mass-to-light 
ratios inferred in WBC09! from which M* and Mjyn are derived. 



2006 WBC09[ ). The comparison of our stellar and dynamical mass 
TFRs to the baryonic TFR is therefore meaningful, and we note that 
the |Stark et al.| ( [200 9l TFR is broadly consistent with our spirals. 
The good agreement gives us confidence in our mass models. 

It is crucial to emphasize that at no point in our analysis do we 
do anything that might systematically affect the SOs in the sample 
differently to the spirals. One slight but unavoidable observational 
bias may be introduced by the presence of dust along the major 
axes. If the slit is placed exactly on the major axis, then absorption 
by dust may make it impossible to recover the full line-of-sight ve- 
locity distributions (LOSVDs). Depending on the optical thickness, 
this may truncate the low velocity wing or, in very dusty discs, the 
high velocity peak of the LOSVDs. When Gaussians are fitted to 
these distribution jChung & Bureau|2004^ , this could therefore in- 
crease or, perhaps more likely, decrease the derived line-of-sight ve- 



However, in cases where the dust absorption was strong, the slit was 
simply shifted from the major axis by ~ 2 arcsec to ensure the full 
LOSVD was sampled, and the objects with strong dust lanes were 
selected to be very slightly away from perfectly edge-on ( [Bureau &] 
[FreemanI 19991 ). 

If the rotational velocity is a strong function of z, the height 
above the equatorial plane in the cylindrical coordinate system of 
the galaxy {R,0,z), this intentional shift could itself bias the ve- 
locities and inferred masses in dustier galaxies (i.e. spirals). In 
a future work (Williams et al., in preparation), we will present 
pseudo-integral-field data (i.e. a sparse velocity fields constructed 
from multiple long-slit positions) for six of the galaxies in the 
present sample. In these data we see no evidence that d v/d |z| ^ 
in the disks of five of the galaxies. At most, there is a gradient of 
— 30kms^' kpc^' in the disk of NGC 7123. This is broadly con- 
sistent with observations of planetary nebulae in the edge-on spi- 



ral NGC 891 (Merrifield et al.[[20I0 i and of extraplanar gas (e.g. 
NGC 891, Heald et al.|2007|[Oosterloo et al.|2007l NGC 4302 and 
NGC 5775, Heald e t al.|2007[ the Milky Way, [Levine et al.|2008^ , 
which find a typical falloff in rotation dv/d|z| ~ — 20kms^' kpc^ 
In the worst case scenario, in which all our spirals (and none of the 
SOs) suffered from absorption strong enough to require shifting the 
slit, the spiral velocities would typically be biased by < — 0.02dex. 
We estimate that the overall effect on our sample of spirals, most of 
which did not require this slit shift, may be to bias their velocities 
by up to ~ — O.OI dex. In practice, the slit was 1.8 arcsec wide and 
rather difficult to align precisely for both spirals and SOs. This im- 
precision, which introduces velocity scatter in both samples, likely 
overwhelms any systematic velocity falloff in the spirals due to in- 
tentionally shifting the slit. 



4 FITTING PROCEDURE 

For each measure of luminosity or mass as a function of x = 
log(vniodel/km s^'), we simultaneously fit a straight line of the 
form y = a{x — 2A) + b to the spiral galaxies and >■ = a(x — 2.4) + 
b + S to the SOs. The two straight lines are therefore constrained to 
have the same slope, a, but they have zero-points which differ by 
5. In separate fits to the two samples, we find that their slopes are 
consistent within the uncertainties. Constraining them to be equal 
therefore significantly simplifies the interpretation of our results in 
terms of zero-point evolution. The x = 2.4 reference value is de- 
fined to minimize the covariance between errors in a and h (e.g. 
[Tremaine et al.|2002| ). Its choice does not affect the results of this 
work, in which were are mainly interested in the zero point offset 
5. 

To find the optimal combination of a, b and 8 for a given scat- 
ter, we define and minimize the figure of merit 

WsaSb 



X=L -2{,y/-[a(x,-2.4)+i]}2 



1=0 "; 

A'so J 



(3) 



^— {yi-[a{xi -2.4) +b + 8]y 



where A'saSb is the number of spirals, A'so the number of SOs and 

(4) 



7 9 9 9 9 

(7/ = CT-,+a^cr-, + CT4- 
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The extra scatter C7;„ 



v^ - 
A- red 



is iteratively adjusted to ensure that 

2 

1. 



X^ 



(5) 



(A'Sa,Sb+A'sO-3) 

If the observational errors for each galaxy C!>_,- and a,./ are well- 
estimated then (Tint is the intrinsic, astrophysical scatter in the TFR. 
We also define a total rms scatter CTtot> essentially a weighted 
mean of the scatter of the data about the fits (and therefore in the 
same units as the y-axis): 

^SaSb 1 

X ^{.y,--Kx,--2.4)+il}2 



Cfnt = 



i=0 



A'so 

+ 1 
1=0 



{yr 



-2.4)+b+d]y 






lf5i/<^f 



(6) 



(7) 



where A'g^j is the total number of both spiral and SO galaxies. Com- 
paring Otot with CT[nt gives an idea of what fraction of the scatter is 
intrinsic and what fraction is due to measurement errors. We mini- 
mize equation pi using the MPFIT package ( |Markwardt|2008^ . 

Following previous analyses of the Tully-Fisher relation, we 
also fit the 'inverse' relation, i.e. we regress the observed rota- 
tional velocities rather than the luminosities or masses onto the 
model (see, e.g., TPOO, [\^rhegenl[200T| fPizagno et al.||2007| and 
references therein). We rewrite the spiral and SO straight lines as 
X = Ay + B + 2.4 and x = Ay + B + 2.4 + A and minimize 

A'SaSb 1 

-{Ayi+B + 2.4)f 



X^ = 



1 



'■=0 £/■ 



1 






where 



and 



^2 I /i2„2 



btot M. , ^ ■ 



(8) 



(9) 



(10) 



To compare these values to those determined using the forward rela- 
tion, we use the fact that a best-fitting inverse relation has an equiv- 
alent forward slope 1/A, zero-point —B/A, offset —A/A and y-axis 
intrinsic and total scatters AQnt and Aijtot- Our implementation of 
the fitting procedure described above is publicly availablaM 

Forward and inverse fitting is only symmetric if there is no in- 
trinsic scatter (aim = Snt = 0) ' Tremaine et al.|2002| l. This is not 
generally the case and the slopes can be very different. In analyses 
of the TFR with observed data, it has been shown that there is a 
significant bias in the slope of the forward line of best fit l |Willick| 
|1994 l, which is why the inverse relationship is usually preferred. 
The figures in this paper use this inverse fitting approach, but to 
aid comparison with future work, we present both forward and in- 
verse fits in Table [2] The slopes of the forward and inverse fits to 
our data are indeed discrepant, often to the extent that they do not 
lie within each other's error bars (especially at 6-band). We are for- 
tunate, however, that the choice of whether to use the forward or 
inverse relationship does not affect our conclusions, which depend 
entirely on the zero-point offset of the SO TFR and the total scatter. 



5 DISCUSSION 

5.1 The TFR zero-point offset at ^s-band and i?-band 

While the TFR zero point and offset are formally sensitive to the 
choice between forward or inverse regression, in practice they are 
almost unchanged by this choice. In the discussion that follows, we 
adopt the means of the forward and inverse values and their errors. 

We detect a statistically significant difference between the 
zero-points of the TFRs of spiral and SO galaxies (see Table|2]l. At 
a given Vmodeb local SOs are 0.53 ± 0. 15 mag fainter at A^^-band and 
0.68 ± 0.22 mag fainter at B-band than local spirals. The systematic 
uncertainty is due to the slit positioning issue described at the end of 
SectionIS] We emphasize again that the B-band offset we measure is 
dubious because of the very uncertain (and very large) internal ex- 
tinction corrections at optical wavelengths in edge-on spirals. The 
A'j-band value is, however, robust. 

We estimate that the intentional misalignment of the slit in the 
dustiest galaxies may introduce a bias of up to ~ —0.01 dex in the 
observed velocities of spirals (see the end of SectionlSJi. In principle, 
this may bias our measurement of the A'^-band offset to be too large 
by ~ 0.07 mag. This is less than uncertainty in the intercept of the 
spiral TFR (iO.lOmag) or the offset of SO TFR (±0.15 mag), so 
we do not discuss it further. 

We restate that our fiducial TFRs use Vmojel as the measure 
of rotation. We emphasize, however, that the offsets we measure 
between the spiral and SO TFRs are completely insensitive to this 
choice. Using i'(ji-ift, there is anoffsetat ATj-bandof 0.45±0.18mag. 
Using Vgas, there is an offset of 0.48 ± 0.24 mag. Unless all three of 
these measures of rotation are flawed then our principal result, the 
measurement of the offset of the SO TFR, is reliable. (We neglect 
vhi because it is based on spatially unresolved data and sometimes 
unphysically discrepant from the other measures of rotation. See 
Section[23|) 

Since high redshift spirals are thought to be the progenitors of 
local SOs, we would really like to compare our local SO TFR to the 
high redshift spiral TFR. As we emphasize in Section [2] however, 
the comparison of TFRs derived using different measures of rota- 
tion is uncertain at best, so the zero-points of high-redshift TFRs 
cannot be compared to our local results in a meaningful way. In ad- 
dition to the uncertainties in comparing the velocities of different 
tracers in the local universe, the move to high-z raises the possibil- 
ity that the intrinsic shapes of rotation curves have changed, which 
we cannot rule out. Provided, however, that these differences have 
been correctly accounted for by other authors, the relative shifts 
between the local and high-z TFRs should be fairly insensitive to 
the measure of rotation used. We therefore make that unavoidable 
assumption in order to proceed. 

[Conselice et al.H2005) , |Flores et al.|p006| > and |Kassin et al.| 



' http://purl.org/mike/mpfitexy 



poo? ) find no evidence for evolution in either the slope or zero- 
point of the A'-band spiral TFR from the local universe to z ~ 1 . 
However, |Puech et al.| ( |2008) detect a dimming of 0.66 mag at Kg- 
band from z = to z ~ 0.6. We first discuss the possibility that the 
spiral TFR has not evolved with redshift. If this is correct, our re- 
sults imply that local SOs are 0.53 ±0.1 5 mag fainter at ATj-band for 
a given rotational velocity than their presumed spiral progenitors. 

To get a feeling for how long such a fading would take, we 
follow [BAM06) by using the |Bruzual & Chariot] pOOS} stellar pop- 
ulation synthesis code to assign an approximate timescale to such 
a fading under various star formation scenarios. We assume solar 
metallicity and a Chabrier initial mass function ( |Chabrier||2003| ) 
and use an updated prescription for the thermally-pulsing asymp- 
totic giant branch ( [Marigo & Girardi|[2007| Chariot 2009, pri- 



© 2010 RAS, MNRAS OQO.fTpT] 



Table 2. Parameters of the best fitting TFRs, wiiere x - 
and y = a{x- 2.4) +b + S for SOs. 
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- log(vniodel/kms ), >' = a(.v — 2.4) + b for spirals 



y 


a 


b 


5 


0\m 


Clot 






Forward regression 








Mk, 


-8.15±0.76 


-24.53±0.11 


0.55±0.15 


0.27 


0.37 


Mb 


-7.67±1.02 


-20.82±0.14 


0.73±0.20 


0.39 


0.50 


log(M*/M(:;) 


3.45±0.26 


11.14±0.04 


-0.22±0.05 


0.06 


0.13 


log(Mdy„/Me) 


3.34±0.27 


11.23±0.04 
Inverse regression 


-0.20±0.05 


0.07 


0.13 


M;,, 


-9.13±0.87 


-24.49±0.12 


0.51±0.16 


0.28 


0.39 


Mb 


-9.81±1.26 


-20.76±0.16 


0.64±0.23 


0.44 


0.56 


log(A/,/Me) 


3.57±0.27 


11.13±0.04 


-O.21±0.O5 


0.06 


0.13 


log(Mdy„/M0) 


3.52±0.28 


11.22±0.04 


-0.19±0.05 


0.08 


0.13 



Notes. (Tii,| is the intrinsic scatter required to yield x^eA = 1.00 ±0.01 and c7|ot is the total scatter (see 
equationpl. a is in units of y/.v, all other quantities are in units of y, i.e. magnitudes or dex of solar masses. 



vate communication). Following a constant star formation episode 
lasting 5 Gyr, the synthetic stellar population takes 0.9^q g Gyr to 
fade by 0.53 ±0.15 mag at Ar^-band and 0.2+[]-^ Gyr to fade by 
0.67 ± 0.21 mag at S-band. If we add an instantaneous 'last gasp' 
burst of star formation comprising 10% of the mass of the galaxy at 
the end of the 5 Gyr episode pedregal et al.|2008^ , the timescales 
increase to 1 .4+^2 Gyr at ATj-band and 0.5+o j Gyr at B-band. 

For plausible star formation histories, the timescales implied 
by the offsets at the two wavelengths are inconsistent. This is per- 
haps not surprising given the uncertainty in the dust corrections at 
S-band. However, the A'^-band offset is far less susceptible to this 
possible systematic error, and this measurement implies an uncom- 
fortably short timescale since the truncation of star formation of 
< 1 .4 Gyr, coiTesponding to a redshift z ~ 0. 1 . As noted byJBAMOS] 
it would be a surprising coincidence if we were living in an era so 
soon after the SOs in our sample transformed from spirals. Obser- 
vationally, this possibility appears to be ruled out by high-redshift 
observations: while [Dressier et aL] ( |1997| l show that the SO fraction 
has risen at the expense of spirals since z ~ 0.5 (5 Gyr ago), |Fasano| 
|et al.| (2000b find that the present relative fraction of SOs and spi- 
rals was largely already in place in clusters at z ~ 0.2 (2.5 Gyr ago) 
and there was no shortage of SOs in groups at z ~ 0.4 (4 Gyr ago, 
[Wilman et al. 20091. In the absence of strong environmental pro- 
cesses, we naively expect the transformation of field spirals to SOs 
to be gradual, and therefore that present day field SOs would need 
to have begun their transformation even earlier. We therefore con- 
clude that a simple scenario of passive fading of spirals into SOs 
is inconsistent with the TFRs of our sample and the evolution of 
the so/spiral fraction with redshift. This interpretation of the off- 
set of the SO TFR is consistent with that of ,Neistein et al.| ( |1999t . 
With limited data they measured a very similar SO TFR offset to 
us, but with much larger uncertainties, and with respect to a line 
width-based TFR (a comparison we argue is flawed in Section[2j. 

We now discuss the possibility that the zero-point of the TFR 
was fainter at earlier times. If the A"-band spiral TFR has evolved 
with redshift and was, for example, 0.66 mag fainter at z ~ 0.6 than 
the local relation (as found by |Puech et al.|2008) , then present-day 
SOs are approximately as luminous at a given rotational velocity as 
spirals were 6 Gyr ago, i.e. there is no evidence for any evolution 
in the velocity-luminosity plane between local SOs and their pre- 
sumed high redshift progenitors. A simple star formation truncation 
scenario is then even harder to reconcile with our results. 

The above analysis is based on the assumption that an SO of 



a given rotational velocity used to have the luminosity of a high- 
redshift spiral of the same rotational velocity. From this assumption 
we have gone on to consider in isolation the evolution of a single 
broadband colour. The parameter space of possible star formation 
histories is, however, highly degenerate and this is far from the ideal 
way to constrain star formation histories. It is nevertheless clear 
from our sample that something other than simple fading beginning 
at the redshifts at which SOs are first observed (and without subse- 
quent star formation), is transforming spirals to SOs. The evidence 
for recent (and even ongoing) star formation in local SOs from UV 
fKaviraj et al.'2007''Jeong et al.'2009 1 and mid-infrared f Temi et aL] 
^2009a b Shapiro et al. 20101 observations is in fact strong. 



5.2 Stellar and dynamical mass 

We also detect a significant difference between the zero points of 
the SO and spiral stellar mass TFRs. The SOs have around 0.2 dex 
less stellar mass for a given rotational velocity than the spirals. Be- 
cause of the relative gas richness of spirals, the offset in the baryonic 
relation (i.e. stellar mass + gas mass, e.g. |McGaugh et al.|2000| ) is 
likely larger and certainly no smaller. Assuming that Vmodel traces 
the enclosed dynamical mass equally well for SOs and spirals, this 
might naively suggest that SOs have less stellar mass per dynamical 
mass within 2-3 Re (the extent of our kinematic data), and therefore 
contain fractionally more dark matter than spirals. This is generally 
consistent with the idea that SOs are older systems that are found in 
denser environments at the bottom of deeper potential wells. If this 
is indeed the case, using dynamical rather than stellar mass should 
eliminate the offset. The dynamical mass TFR is essentially a plot 
of a one measure of dynamical mass (from a mass model) against 
another (the observed rotational velocity), and since gravity applies 
equally to spirals and SOs there should be no difference. 

Having said that, the difference between the spiral and SO zero 
points persists almost unchanged in the dynamical mass TFR (see 
Table[2l. Our dynamical modelling approach does necessarily make 
approximations and assumptions, but it is unlikely that the models 
could be systematically biased by morphological type. Our stellar 
and dynamical mass spiral TFRs are also consistent with the reli- 
ably calibrated baryonic TFR of |Stark et al.| (|2009[>. This raises the 
possibility that SOs are systematically smaller or more concentrated 
than spirals. Recall that M^y^ = v^R, where R is some characteristic 
radius. If /? is systematically smaller in SOs at a given rotational ve- 
locity, or if the dynamical mass is more concentrated such that more 
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mass is enclosed within a given characteristic radius, then that could 
explain the offset in the dynamical mass TFR. In other words, spi- 
rals and SOs do not form a homologous family. In this scenario, the 
offset of the SO dynamical mass TFR should not be thought of as 
a -0.18 dex offset in mass (as presented in Table |2l, but, because 
of the slope of the TFR, a +0.05 dex offset in velocity. This could 
be due either to SOs being approximately 80% of the size of spirals 
of the same rotational velocity, or their dynamical mass being dis- 
tributed more compactly such that the characteristic radius probed 
by measures like v,nodel and vjrift contains approximately 25% more 
dynamical mass. Again, since SOs are thought to be older, more 
evolved systems with bulge-to-disc ratios that are larger on average 
than those of spirals, this is perhaps consistent with naive expecta- 
tions, but is it true for our sample? In the spirit of Courteau & Rix| 
( |I999[ > and |Courteau et alT]j2007| >, the possibility that the luminous 
component of the dynamical mass is smaller or more concentrated 
in SOs at a given velocity can be tested directly by measuring the 
size of our galaxies. 

We define R^ as described in Section [2^ (where it is used as 
a parameter of the asymmetric drift correction) and R^ as the semi- 
major axis of the ellipse containing half the light at A'^-band (see 
section 4.2 of |WBC09l for a discussion of alternative definitions). 
In an edge-on galaxy with a boxy bulge, such as those in our sam- 
ple, the radial surface brightness profile often contains a plateau or 
even a secondary maximum ( [Bureau et al.|2006) . This means the 
disk scale length is not well defined. Measurement of R^ depends 
sensitively on the ellipticity of the aperture used |WBC09 |>. More- 
over, our data probe a relatively small range in Vmodel- Together, 
these difficulties mean that the constraints we can place on the pa- 
rameters of the size-velocity relation are very weak. We recover a 
correlation i?(j txR^x I'^^^^gj, where y~ 1.5 ±0.5. Significantly, we 
also detect some evidence of a systematic difference between the 
zero points of the size-velocity relations of SOs and spirals. Using 
the simultaneous line fitting approach described in Section H] we 
find that SOs are on average 0. 15 ± 0. 10 dex smaller than spirals at 
a given v,nodel- This result is of weak statistical significance, but is 
consistent with our interpretation. We see no significant difference 
between the mean compactness of the light distribution, C%2^ of the 
SOs and spirals. Cf,2 = ^80 /^e, where R^q is the semi-major axis of 
the ellipse containing 80% of the light. 

The large sample of [Courteau et al. | j2007[ l does not include 
SOs, but extends from Sa to Sc and is made up of galaxies that are 
more favourably inclined for measurements of size. Recalling that 
^Lj-band luminosity is a good proxy for stellar mass (e.g. [Bell &| 
[de Jong|2001[ l, which is a reasonable proxy for dynamical mass in 
the optical parts of disc galaxies (e.g. [van Albada & Sancisi 1986; 
Englmaier & Gerhard|1999[|Palunas &WilliamS|2000,|Weiner et al. 



2001) , [Courteau et al.[ ( [2007 > find a morphological trend in the sense 
required to explain our result: earlier-type spirals are smaller for 
a given A'^-band luminosity than later-types (see their table 3 and 
figure II). At a characteristic A^^-band luminosity of 10 Lq, they 
find that Sas are around 80% the size of Sbs, which are in turn 
around 80% the size of Scs. If this trend extends to SOs, then its 
sense and magnitude are consistent with the offset we see in our 
dynamical mass TFR. 

There is also the possibility that the dark halos (rather than the 
baryonic mass distribution) of SOs are smaller or more concentrated 
than those of spirals. Unfortunately we are unable to quantify this 
with our kinematic data, which are not radially extended enough 
to break the well known degeneracy between halo mass and con- 
centration l [van Albada et al.][1985) . However, the inability of the 
dynamical models presented in WBC09|to quantify the amount of 



halo contraction has no effect on the results presented here. This is 
because, even if one assumes different dark halo shapes, the param- 
eters of the best fitting models always conspire to produce an almost 
unchanging v^ profile in the region constrained by the data. This is 
ultimately the reason why the degeneracy exists in the first place. 
A quantitative illustration of the robustness of the circular velocity 
profiles recovered from dynamical models to the degeneracy in the 
halo concentration is given in, e.g., fig. 9 of [Thomas et al.[(2009^ . 

The interpretation of the systematic offset of the SO TFR to 
higher velocities due to a more concentrated mass distribution or 
smaller size than spirals is consistent with the discussion of the 
shapes of very extended gas rotation curves in |Noordermeer et al.| 
( [2007[ > and [Noordermeer & VerheijenH200 7"'). Their kinematic data 
extend well beyond the optical discs of the galaxies, and allow them 
to study the asymptotic behaviour of the rotation curves. They find 
that the earlier galaxies in their sample of SO-Sab galaxies have ro- 
tation curves that reach a maximum that is greater than their asymp- 
totic velocity and occurs at smaller radii. They argue that this be- 
haviour can be understood in terms of the larger bulge-to-disc ra- 
tios of the earlier-type galaxies in their sample. The scatter of their 
earlier-type galaxies away from the TFR of the late type galaxies 
is reduced by using the asymptotic velocity, which again suggests 
the Tully-Fisher relation is a manifestation of a close connection 
between galaxy luminosity (or mass) and halo mass. 

If the size/concentration argument discussed above is indeed 
the explanation for the difference between the zero points of the SO 
and spiral TFRs as a function of dynamical mass, then the offsets we 
observe as a function of luminosity (and discuss extensively in Sec- 
tion [5]T] in the context of fading relative to the high-redshift spiral 
TFR) are in fact due to differences in the distributions of dynamical 
mass in spirals and SOs. The offset between spirals and SOs in the 
Tully-Fisher (luminosity-velocity) relation may therefore not be a 
property of their stellar populations as argued in BAM06 , but rather 
the result of more fundamental offset between spirals and SOs in the 
mass-size (or luminosity-size) projections of the luminosity-size- 
velocity plane occupied by disc galaxies. Models of SO formation 
in which SOs are end points of spiral evolution should consider this 
possibility. 



5.3 Scatter of the TFRs 

The B-band luminosities we use, which are drawn from HYPER- 
LEDA, are less accurate than our jPj-band magnitudes, which are 
based on a relatively recent, homogeneous set of observations. The 
extinction corrections at S-band are also particularly large and un- 
certain for edge-on galaxies like ours. These uncertainties are prob- 
ably not reflected properly in our error estimates, and this may con- 
tribute to the fact that the scatter of the TFR at S-band (0.53 mag, of 
which 0.42 mag is intrinsic rather than observational) is larger than 
that at A'5-band (0.38 mag, of which 0.28 mag is intrinsic). 

These scatters are broadly consistent with dedicated studies of 
the spiral TFR scatter (e.g. [Pizagno et al.|2007) , which is perhaps 
somewhat surprising given the expected poor quality of our distance 
estimates and the importance and uncertainty of internal extinction 
in edge-on spirals. It should be borne in mind, however, that such 
studies fit a single line to a range of morphological types, while 
we separately fit two lines to two sub-samples, which reduces the 
measured total and intrinsic scatters. If a single line is fitted to the 
composite sample of spirals and SOs then the scatter increases to 
0.50 mag, of which 0.41 mag is intrinsic. 

Because our sample is drawn from a wide range of environ- 
ments, the SOs may have been transformed from spirals via mul- 
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dple channels (or one channel which began at different times in 
different galaxies). If so, then the extent to which SOs have left the 
spiral TFR should vary, and their scatter about the TFR should be 
larger than that of spirals. Our fitting method implicitly assumes 
that spirals and SOs have the same intrinsic scatter. However, our 
analysis also shows that there is no evidence of any difference be- 
tween the intrinsic or total scatter of spiral galaxies about their line 
of best fit and that of SOs about theirs. The absence of increased 
scatter in the SO TFR is problematic for all explanations of the off- 
set (passive fading, environmental stripping, minor mergers and the 
possible non-homology described in Section [F!2| . 

There is no significant difference between the total scatters of 
the stellar mass TFRs (0.13 dex, i.e. 35 per cent) and the scatter at 
ATj-band (0.32 mag, i.e 34 per cent in luminosity). This is proba- 
bly due to the relative constancy of mass-to-light ratios in the near- 
infrared, from which the masses are derived. The difficulty of as- 
signing errors to the total masses of the dynamical models makes 
drawing strong conclusions from the intrinsic scatter of the stellar 
mass or dynamical mass TFRs difficult. 

5.4 Classification of edge-on disc galaxies 

The ability to reproducibly and correctly distinguish between spiral 
and SO galaxies is crucial for this work. Although spirals arms are 
not visible in edge-on disc galaxies, the classification of edge-on 
galaxies as spirals or SOs is clearly specified by the presence of ex- 
tended dust (see, e.g., |Hubble|1936[|Sandage|1961[[de Vaucouleurs| 
|et al.|199 Tl|. It is clear, therefore, that the classification of edge-on 
galaxies is objective and reproducible. 

Although the classification of edge-on disc galaxies is repro- 
ducible, would a galaxy be classified as the same type if viewed 
face-on? To answer this concern, we note that if the effects of dust 
are visible in a face-on spiral then, due to line-of-sight projection, 
they would be even clearer if the same galaxy were viewed edge- 
on. For the simple distinction between spiral and SO galaxies (rather 
than more detailed division into the Sa-Sd classes based on bulge 
size and the tightness of spiral arms), edge-on orientation is thus 
arguably the optimal viewing angle. This is also true of the dis- 
tinction between SOs and ellipticals, where a face-on orientation 
makes it extremely difficult to detect the featureless stellar disc of 
an SO. Having said that, it is clear that if an edge-on galaxy has just 
enough dust to be classified as a spiral, it would probably (but per- 
haps wrongly) be classified as an SO if viewed face-on. This will 
occur only for limiting cases, however, and it is clear that the ma- 
jority of edge-on classifications are correct. 

Having established that a meaningful morphological classifi- 
cation of edge-on disc galaxies is possible, we now ask whether this 
simple spiral-SO division based on the presence of dust is dynam- 
ically significant. To test this, we randomly divided the complete 
sample of 28 galaxies into two sub-samples of 14 galaxies each and 
then measured the offset 8 between the zero points of the TFRs of 
the two random sub-samples. We repeated this procedure 50,000 
times to build up an idea of how the offset between randomly se- 
lected sub-samples is distributed, and how exceptional the adopted 
morphological classification is. This idea is a simplification of the 
non-parametric test presented in the context of comparing local and 
high-redshift TFRs by K oen & Lombard] (2009). Our results are 
presented in Fig.B] It is clear that the offsets observed between our 
SOs and spirals cannot to be due to chance. In all cases (y = M^^, 
Mb, Mi, or M^y^, and for both forward and inverse regressions), the 
observed offset is at least three standard deviations away from the 
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Figure 5. Histograms showing the TFR zero-point offset 8 between two 
sub-samples drawn randomly (with no regard to classification) from our full 
sample of 28 galaxies. The inverse regression fitting procedure was repeated 
50,000 times to construct the histograms. From top to bottom, we show the 
offsets for the X^-band, 6-band, stellar mass and dynamical mass TFRs as a 
function of Vmodd- The vertical lines show the value of the zero-point offset 
measured for the spiral and SO sub-samples used here. 



mean of the offsets between randomly selected sub-samples. This 
result clearly demonstrates that a purely morphological classifica- 
tion of edge-on early-type disc galaxies, based only on the presence 
of dust, divides them into two truly kinematically distinct classes. 



5.5 Comparison to previous work 

The absolute location of our SO TFR is almost identical to that mea- 
sured in BAM06 at both Kg and 6-band. They used Vdijft rather than 
Vmodeh SO this is a reflection of the consistency of these two mea- 
sures of rotation. If |B AMOS] could have used our spiral TFR as a 
reference zero point, they would have measured the same offset as 
us. However, they used the TFRs of TPOO. As a result, our offsets 
are in the same sense (SOs are fainter than spirals) but significantly 
smaller in size than those measured in[BAM06](l .0±0.4 mag at Ks, 
1 .6±0.4 mag at fi-band, where we have removed their correction of 
the TPOO zero point). Because the zero-point offsets we measure are 
smaller than those in[BAM06| the timescale for a simple synthetic 
stellar population with a plausible star formation history to fade is 
shorter too. There are at least two reasons for the larger offsets mea- 
sured in' BAMO^ which could combine to account for the entirety 
of the difference, making our results consistent. 

The first is the possible intrinsic difference between the zero 
point of the spiral TFR we compare our SOs to in the present work, 
and the spiral TFR used in jBAMO^ Our spiral TFR is constructed 
using Sa and Sb galaxies, while that of TPOO, which is used by 
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|BAM06[ is calibrated using later types (mostly Sbc and later). As 
shown by |Masters et al.|J2008^ , in the velocity domain of the present 
work ((v) ~ 250 kms^'), Sa spirals are ~ 0.4 mag fainter than Sb 
and Sc spirals of the same rotational velocity. On the other hand, 
[Courteau et al.| ( (2007| l estimate that Sas are just 0. 1 mag fainter than 
Sbs. Indeed, with our much smaller statistics, we see no evidence of 
a zero point variation between the Sa and Sb galaxies in our sample. 

The second possible effect is the systematic bias introduced 
by comparing an SO TFR derived from stellar kinematics to a spiral 
TFR derived from global H I line widths. Comparing the methods 
shows that Sa-Sb galaxies are measured to be rotating ~ 0.08 dex 
slower in H I than in their stellar kinematics (see Fig.[2](. At a typ- 
ical ^5-band TFR slope of ~ 8, this corresponds to a magnitude 
shift of ~ 0.6 mag. However, as we emphasized in Section [2] ap- 
plying a constant shift to correct for the difference between stellar 
and H I rotation measures is dubious because of the large scatter in 
this difference (admittedly less so for spirals). 

It is easy to see how morphological differences between galax- 
ies in the reference spiral samples could combine with the system- 
atic bias introduced by the use of different measures of rotation to 
account for the discrepancy between the fiducial offset presented 
here (0.53 ± 0. 15 mag) and that in BAM06 (1 .0 ± 0.4 mag). 

Another possible source of bias in the offset measured in 
|BAM06| is the fact that the TPOO spiral TFR is constrained by 
observations of a sample o f relatively low mass spirals with v ~ 



150kms-l, while |BAM06 



i[s SOs have a mean v ~ 250 km s^. If 
the slope of the TFR changes with mass, as has been suggested (e.g. 
[Peletier & Willner 1993 ; Noordermeer & Verheijen 2007 ), then the 
fact that the SOs and spirals in our sample lie in the same velocity 
domain is a crucial advantage. 

In any case, because the rotation measure bias is unphysical 
and unrelated to the formation and evolution of SOs, we argue that 
the method we have employed here should be preferred if the goal 
is to determine how much fainter local SOs are than local spirals. 

Finally, we note that |BAM06| found a much larger scatter in 
their SO TFR than we do here, ~ 0.9 mag at A^j-band. This may 
be explained by the relatively homogeneous nature of our data 
compared to the multiple sources from which [BAM06| drew theirs, 
and perhaps to the larger velocity and luminosity ranges their data 
probe. 



6 CONCLUSIONS 

We have demonstrated that the comparison of TuUy-Fisher relations 
derived from H I line widths, ionized gas PVDs, stellar kinematics 
corrected for asymmetric drift or the circular velocity profiles of 
dynamical models is influenced by systematic and uncertain biases 
introduced by the different measures of rotation used. We therefore 
argue that to constrain the relative locations of the spiral and SO 
TFRs, the same tracer and methods must be used for both samples. 
In practice, because of the paucity of extended, undisturbed H I and 
ionised gas in SOs, this means one must use stellar kinematics or 
dynamical models. 

In this work we used the circular velocity profiles of mass 
models to construct TFRs for 14 Sa-Sb spirals and 14 SOs, elim- 
inating the biases introduced by mixing measures of rotation. The 
circular velocity curves are those of models derived by solving the 
Jeans equations for mass models comprising an axisymmetric stel- 
lar component and a spherical NFW halo ( |WBC09| l. The parame- 
ters of the models are constrained by observed long-slit major-axis 
stellar kinematics l|Chung & Bureau]|2004|>. We characterized the 



circular velocity profile for each galaxy by single numbers, v^Qdei, 
by taking the average in the flat part of the observed rotation curve. 

By simultaneously fitting two offset relations with a common 
slope to this spirals and SOs, we find that SOs are systematically 
0.53 ± 0. 15 mag fainter at A'^-band than local Sa-Sb spirals of the 
same rotational velocity. This measurement is almost identical if 
we use estimates of the rotational velocity derived from ionized gas 
PVDs or stellar kinematics corrected for asymmetric drift. 

If the high-redshift spiral TFR has the same zero point as the 
local spiral TFR, this is inconsistent with the observed evolution of 
the spiral/so fraction with redshift and a simple scenario in which 
star formation in the spiral progenitors of SOs was truncated at some 
time in the past. More complex star formation histories or even on- 
going star formation in SOs may be the explanation. An alternative 
interpretation is revealed by the stellar mass and dynamical mass 
TFRs. The TFR offset persists as a function of both stellar and dy- 
namical mass, and we show that this may be evidence of a small 
(10-20%) but systematic contraction of spirals as they transform to 
SOs. This is consistent with the trend with morphological type of the 
size-luminosity relation in the local universe jCourteau et al.|2007) . 
If, on the other hand, the zero point of the TFR has dimmed from the 
present day to high redshift by ~ 0.5 mag, then the putative trans- 
formation from spirals to SOs involves essentially no movement in 
the velocity-luminosity plane. 

It seems clear that SOs are not primeval objects, but are an end 
point of spiral evolution. The processes responsible for this evolu- 
tion can perhaps be accelerated by the environment. The variation 
of the zero point of the TFR with galaxy type and other parame- 
ters is just one approach among a number available that will allow 
us to constrain SO formation and evolution. It is uniquely powerful 
because it encodes information about the amount and distribution 
of dynamical mass of a galaxy. Modelling work complementary to 
the observations presented here is underway ( |Trujillo-Gomez et al.| 
|20I0[ |Tonini et al.||2010| >. One should be sure that the inevitable 
compromises introduced by high redshift observations do not intro- 
duce biases similar to those ones discussed in Section|2] We should 
also not forget the crucial role that spectroscopy has to play, espe- 
cially in constraining the star formation history. 
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APPENDIX A: DEFINITIONS OF DYNAMICAL MASS IN 
DISC GALAXIES 

In Section [3] we discussed two ways to define or characterize the 
dynamical mass of a galaxy. The first is to multiply the /f^-band 
luminosity L by the dynamical mass-to-light ratio (M/L)/;-^ nohalo 
(taken from |WBC09| for our galaxies). As we noted, this method is 
relatively straightforward, and we adopt it as our fiducial method in 
the main body of this paper (denoting it Mdy,,). This quantity does 
not have a precise physical meaning, but here we show it is equal, 
to a good approximation, to IM^n, where Mj « is the dynamical 
mass contained within the three-dimensional half-light volume. 

We estimate M^n by inverting the expression for the circular 
velocity of the mass model at the rj «. Ror a spherical system (i.e. 
flattening q = c/a =1) this gives 

■■?(. 



M 



111 = ni2^ 



'^i/i 



)/G, 



(Al) 



where vdrin) is the circular velocity at r; m. Note that rjM is not 
the projected half-light (effective) radius, R^. Ror a wide range of 
radial profiles, rj/j ~ l33Re l [Hemquist| 1990} |Ciotti|[T99T| [Wolf 
|etal.|2010t . 

Disc galaxies are not spherically symmetric, so a more use- 
ful definition in these cases is the dynamical mass within an oblate 
spheroid of semi-major axis rj m, which we denote Mj m q- For ^ 
given circular velocity, M^n > M^n q- This is because the mass 
that must be spherically distributed to maintain a given circular ve- 
locity is larger than that necessary to achieve the same circular ve- 
locity in the midplane of a flattened system. The ratio Mj m/Mi m « 
is simultaneously a function of the radial profile of the mass dis- 
tribution and q. Ror a full analytic discussion see, for example, 
the comparison between the circular velocities of oblate spheroids 
as a function of q and mass profile in section 2.5.2 of Binn ey &] 
[Tremaine] j2008| ). For our purposes, an approximate estimate of 
^Xlll^lll q is sufficient. Using a range of analytic potentials, we 



examine the behaviour of the ratio Mi 



/2/M1/. 



2,? 



as a function of 



q in Fig. |A1| The circular velocity of each flattened model is con- 
veniently calculated numerically by flattening a multi-Gaussian fit 
to the mass distribution and using standard integrals to determine 
the circular velocity of the arbitrarily flattened distribution jCap-| 
|pellari|2002|>. Disc galaxies typically have q ~ 0.2 and potentials 
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Figure Al. M\j2lM^j2,q for a selection of potentials as a function of the 
flattening of the potential q. The behaviour of the |Hemquist| jT990l pro- 
file is shown in solid purple. The other lines are Sersic profiles of order 
n=\ (exponential), 2.5 (typical classical bulge-pseudobulge transition), 4 
(de Vaucouleurs) and 10. The de Vaucouleurs and Hernquist lines are well 
described by the fitting formula Mi/2/Mi/2,^ = 1.38-0.52^ + 0.14^^. 



of the form exp(— iJ/i?e)''"> where 1 < n < 4. In such galaxies, 



1.3 < A/j rt/Mj M n ^ 1-5, and so from equation I Al i 



M 



l/2,^~''l/2V?(''l/2)/l-4G. 



(A2) 



Furthermore, rj « 



1 .33i?e> where in non-spherical systems the ob- 
servable R^ is defined as the semi-major axis of the projected half- 
light ellipse. For disc galaxies we can therefore write 

,2/ 



M 



\ll.c 



iRevi{l33Re)/G. 



(A3) 



The circular velocity at 1.33 Re can easily be estimated to an 
accuracy consistent with the derivation of this formula from dynam- 
ical modelling, stellar kinematics corrected for asymmetric drift or 
resolved gas kinematic observations. For typical values of q and 
the analytic potentials explored here, this expression is accurate to 
about 15 per cent. Of course in real disc galaxies, the flattening q 
is a function of radius, the potential is not described by a single 
Sersic or Hemquist profile and even R^ is not simple to measure, so 
equation iJASl is probably a little less reliable in practice. 

Nevertheless, when doubled, Mj m q closely matches fiducial 
estimates of the dynamical mass of our sample galaxies, i.e. the 
product of their luminosities L and dynamical mass-to-light ratios 
(M/L)nohalo5 inferred in this case from stellar kinematic data ex- 
tending to 2-3 Rg and A'^-band photometry in |WBC09| We show 
this comparison in Fig. |A2| We use v^iodel. the circular velocity of 
the |WBC09| mass models in the flat region of the observed rota- 
tion curve, as v^ in equation l|A3j, but in cases where a reliable 
asymmetric drift coiTection is possible at 1.33 R^, v^rift gives almost 
identical results. The excellent agreement between 2Mj m « and 
L X (M/L)„ohalo m''y be a coincidence. The two measures would 
be systematically offset (and perhaps more scattered) if our kine- 
matic data extended significantly further in radius (and the dynam- 
ical mass-to-light ratios we measured were therefore larger). Cru- 
cially, however. Fig. | A2| demonstrates that, for the purposes of the 
present work comparing the TFRs of spirals and SOs, the choice of 



Figure A2. A comparison of the two estimates of the dynamical mass for 
the 28 galaxies in the present sample. The horizontal axis is 2Mi m.o as de- 
fined in equation JA3j, incorporating an approximate correction for flatten- 
ing. The vertical axis is the A'j-band luminosity L multiplied by the dynami- 
cal mass-to-light ratio {M / L) ^^ ^„t^^i„ as presented in WBC09 The symbols 
are as in Fig.l2] The solid line is the line y = x. The two estimates agree on 
average to better than 2 per cent with 13 per cent rms scatter and no evidence 
of a systematic difference between spirals and SOs. 



definition of dynamical mass does not affect spirals and SOs differ- 
ently. 
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